Hypernuclei, dibaryon and antinuclei production in high energy heavy ion 
collisions: Thermal production vs. Coalescence 1 ^ 
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Abstract 

We study the production of (hyper-)nuclei and di-baryons in most central heavy Ion collisions at energies of 
Eiab = 1 - 160A GeV. In particular we are interested in clusters produced from the hot and dense fireball. The 
formation rate of strange and non-strange clusters is estimated by assuming thermal production from the intermediate 
phase of the UrQMD-hydro hybrid model and alternatively by the coalescence mechanism from a hadronic cascade 
model. Both model types are compared in detail. For most energies we find that both approaches agree in their pre- 
dictions for the yields of the clusters. Only for very low beam energies, and for di-baryons including E's, we observe 
considerable differences. We also study the production of anti-matter clusters up to top RHIC energies and show 
that the observation of anti- 4 //e and even anti- 4 He is feasible. We have found a considerable qualitative difference 
in the energy dependence of the strangeness population factor Rh when comparing the thermal production with the 
coalescence results. 
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1. Introduction 

Relativistic heavy ion collisions are an abundant 
source of strangeness. As strange quarks have to be 
newly produced during the hot and dense stage of the 
collision, they are thought of carrying information 
on the properties of the matter that was created fljj] . 
Together with other probes like the elliptic flow and 
jet quenching, the enhancement of strange particle 
production is discussedjl SJ1J1 0, 0, S @,M M El , 
Ll! Q OS H HE EI (HElllll] as a possible signal 
for the creation of a deconfined phase. 
Although abundantly produced, the strong interactions 
of strange hadrons are not well understood. Such 
interactions are not only important for the description 
of the hadronic phase of a heavy ion collision but 
also play an important role for the description of dense 
hadronic matter. In this context hyperon interactions are 
key to understand the phase structure of QCD at large 
densities and the interior of compact stars. One way to 
tackle the problem of hyperon interactions is to study 



the formation of hyperclusters and/or hypernuclei. 
Hypernuclear physics offers a direct experimental way 
to study hyperon-nucleon (YN) and hyperon-hyperon 
(YY) interactions (Y = A, 2, S, Q). The nucleus serves 
as a laboratory offering the unique opportunity to study 
basic properties of hyperons and their interactions. 
Even the confirmation or exclusion of the existence for 
such objects can be used as an input for models that try 
to describe hyperonic interactions. 
More exotic forms of deeply bound objects with 
strangeness have been proposed [22] as states of matter, 
either consisting of baryons or quarks. The H di-baryon 
was predicted by Jaffe 12311 and later, many more 
bound di-baryon states with strangeness were proposed 
using quark potentials 124 12511 or the Skyrme model 
B26I1 . However, the non-observation of multi-quark 
bags, e.g. strangelets is still one of the open problems 
of intermediate and high energy physics. Lattice 
calculations sug gest that the H-dibaryon is a weakly 
unbound system [27], while recent lattice studies report 
that there could be strange di-baryon systems including 
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S's that can be bound B28I1 . Because of the size of these 
clusters lattice studies are usually very demanding 
on computational resources and have large lattice 
artifacts, it is not clear if Lattice QCD predicts a loosely 
bound H-dibaryon or if it is unbound B29l l3Cl l3l[|32ll . 
An experimental confirmation of such a state would 
therefore be an enormous advance in the understanding 
of the hyperon interaction. 

For completeness we also include in our analysis a 
hypothetical NA di-baryon with mass 2.054 GeV (see 
Table 1), a weakly bound state of a A-hyperon and 
a neutron. The search for such an exotic object is 
underway at GSlQ. 

Hypernuclei are known to exist and be produced 
in heav y I on collisions already for a long time 
1 34, I35J \36L I37I1 . The recent discoveries of the first 
anti-hypertriton [38] and anti-a 1 39] (the largest anti- 
particle cluster ever reported) has fueled the interest 
in the field of hypernuclear physics. Metastable exotic 
multi-hypernuclear objects (MEMOs) as well as purely 
hyperonic systems of A's and H's were introduced in 
1 40J |4l|] as the hadronic counterparts to multi-strange 
quark bags [42, 43]. 



Hypernuclear clusters can be produced and studied 
in various experimental setups, e.g. from proton or anti- 
proton induced reactions [440 as well as pion and kaon 
beams JH, 0, H3, H, l49tl. In this work we will focus 



on the production of hypernuclei in high energy colli- 
sions of Au+Au ions |50]. In such systems strangeness 
is produced abundantly and is likely to form clusters 
of different sizes. Our aim is to determine which pro- 
cesses are most efficient in searching for hypernuclei 
including exotic ones. Presently, we can discriminate 
two distinct mechanisms for hypercluster formation in 
heavy ion collisions. First, the absorption of hyperons 
in the spectator fragments of non central heavy ion col- 
lisions. In this scenario one is interested in hyperons 
which propagate with velocities close to the initial ve- 
locities of the nuclei, i.e., in the vicinity of nuclear spec- 
tators I51IJ52I.I53I.I54I1 . The hyper-systems obtained here 
are rather large and moderately excited, decaying into 
hyperfragments later on 1154 15511 . Alternatively, (hyper- 
)nuclear clusters can emerge from the hot and dense 
fireball region of the reaction. In this scenario the clus- 
ter is formed at, or shortly after, the (chemical-)freeze 
out of the system. A general assumption is, that these 
clusters are then formed through coalescence of differ- 
ent newly produced hadrons B56I1 . To estimate the pro- 
duction yield we can employ two distinct approaches 
which allow us to estimate the theoretical uncertainties 
associated with different treatment of the process. First 
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Table 1 : Properties of all considered multibaryonic states 

we use a hadronic transport model to provide us with 
the phase space information of all hadrons produced in 
a heavy ion collision. This information then serves as 
an input for a coalescence prescription. On the other 
hand it has been shown 11571 I58LI59I1 that thermal models 
consistently describe the production yields of hadrons 
(and nuclei 16011 ) very well. We can therefore assume 
thermal production of clusters from a fluid dynamical 
description to heavy ion collisions. 
Both approaches differ significantly in their assump- 
tions and one would expect to obtain different results, 
depending on the method used. Hence it has been pro- 
posed (e.g. see [61, 62]) that the yield of an exotic 
hadronic state may depend strongly on its structure. The 
purpose of this paper is therefore to, comprehensively, 
compare hypernuclei and di-baryon production from a 
coalescence and thermal/hydrodynamical approach, and 
interpret the differences. One particular important point 
is that we deliberately compared two distinctively differ- 
ent models to explore the robustness of our predictions. 
In this way we can estimate systematic differences in- 
troduced by the two models features, for example a dif- 
ference in the baryon stopping or hyperon phase space 
distributions. 



2. Thermal production from the UrQMD hybrid 
model 

The hybrid approach used in this work is based 
on the integration of a hydrodynamic evolution into 



the UrQMD transport model [63], |6J, |65J- During 
the first phase of the evolution the particles are de- 
scribed by UrQMD as a string/hadronic cascade. Once 
the two colliding nuclei have passed through each 
other the hydrodynamic evolution starts at the time 



tstart - 2R/ y/yc.,,,. — 1, where y cm . denotes the Lorentz 
factor of the colliding nuclei in their center of mass 
frame. While the spectators continue to propagate in 
the cascade, all other particles, i.e. their baryon charge 
densities and energy-momentum densities, are mapped 
to the hydrodynamic grid. By doing so one explicitly 
forces the system into a local thermal equilibrium 
for each cell. In the hydrodynamic part we solve the 
conservation equations for energy and momentum 
as well as the net baryon number current, while for 
the net strange number we assume it to be conserved 
and equal to zero locally. Solving only the equations 
for the net baryon number is commonly accepted in 
hydrodynamical models, although we have shown 
in earlier [66] publications that net strangeness may 
fluctuate locally. It is planned to also implement an 
explicit propagation for the net strange density. 
The hydrodynamic evolution is performed using the 
SHASTA algorithm [67]. At the end of the hydrody- 
namic phase the fields are mapped to particle degrees 
of freedom using the Cooper- Frye equation [68] with 
the properties of the clusters, which serve as input for 
the computation, being listed in Table [1] The transition 
from the hydrodynamic prescription to the transport 
simulation is done gradually in transverse slices of 
thickness 0.2 fm, once all cells in a given slice have 
an energy density lower than five times the ground 
state energy density (see also |693). The temperature 
at hb - which corresponds to such a switching 
density is roughly T — 170 MeV which is close to what 
is expected to be the critical temperature. Detailed 
information of the transition curve in the phase diagram 
can be found in 16311 . In this work we neglected final 
state interactions of the clusters produced. This can 
be justified, as previous works have shown that final 
state interactions reduce e.g. the deuteron yield by only 
about 20% H. 

For an extensive description of the model the reader is 
referred to 11631 16911 . 



3. Coalescence from the DCM-QGSM approach 

Another model used to describe the dynamical stage 
of the reaction is the intra-nuclear cascade model de- 

(We refer to it as the 



veloped in Dubna 17 lL 17211 . 
Dubna Cascade Model - DCM.) The DCM is based on 
the Monte-Carlo solution of a set of the Boltzmann- 
Uehling-Uhlenbeck relativistic kinetic equations with 
the collision terms, including cascade-cascade interac- 
tions. For particle energies below 1 GeV it is suffi- 
cient to consider only nucleons, pions and deltas. The 
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Figure 1 : Mass dependence of calculated invariant yields of light frag- 
ments and hyperfragments produced in central Au+Au collisions at 
1 1.5 A GeV/c compared with experimental data [73] for Au + Pb col- 
lisions. The lines are empirical interpolations of the results. 



model includes a proper description of pion and baryon 
dynamics for particle production and absorption pro- 
cesses. At energies higher than about 10 GeV, the 
Quark-Gluon String Model (QGSM) is used to describe 
elementary hadron collisions. The QGSM considers the 
two lowest SU(3) multiplets in mesonic, baryonic and 
antibaryonic sectors, so interactions between almost 70 
hadron species are treated on the same footing. The 
above noted two energy extremes were bridged by the 
QGSM extension downward in the beam energy II72I1 . 

For the present study the coalescence model has been 
modified in comparison with its initial formulation in 
117 811 . As usual, the coalescence model forms a deuteron 
from a proton and a neutron produced after the cascade 
stage of reaction if their relative momenta are within a 
sphere of radius pc, comparable to the deuteron's mo- 
mentum. The same momentum criterion can be used 
to describe formation of tritons, 3 He, and a-particles. 
In particular, the parameters pc(d)-9Q , /?c(0-108 , 
/?c( 3 He)=108 , and p c (a)-H5 (MeV/c) were adopted 
to reproduce the experimental data [78]. An approach 
disregarding the spacial coordinates of nucleons can be 
justified only for collisions with moderate energy depo- 
sition in nuclei since the region for final state interaction 
is small enough. However, this is not the case for cen- 



tral heavy ion collisions. Here we assume that the coa- 
lescence criterion used to form the composite particles 
includes the proximity of nucleons both in the momen- 
tum and coordinate space. The coordinate coalescence 
parameters are determined by the relation re - fr/pc, 
with the same values of pc as were used in [ 78] . As a 
first approximation we use the same coalescence param- 
eters for both conventional fragments and hyperfrag- 
ments. An example of the calculated invariant yields 
of the fragments produced in the central Au + Au col- 
lisions at projectile momentum 11. 5 A GeV is shown 
in Fig. [1] One can understand that at this energy the 
coalescence model reproduces qualitatively the experi- 
mental data for conventional fragments. The fragments 
yields fit very close to exponential dependence with a 
penalty factor of approximately 50 for each nucleon 
added in agreement with the data. Due to the fact that 
the same coalescence parameters were used a similar 
penalty factor is obtained for hyperfragments, which is 
supplemented by additional suppression if the neutron 
is replaced by a A. 

For the following results we fixed the coalescence pa- 
rameters as described, with a fit to the data at 11. 5A 
GeV, and assume that they do not change with beam en- 
ergy. This allows us to predict cluster production over a 
wide range of experimental setups. 



4. Results 
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Figure 2: Yields per event of different di-baryons in the mid rapidity 
region (|y| < 0.5) of most central collisions of Pb+Pb/Au+Au. Shown 
are the results from the thermal production in the UrQMD hybrid 
model (lines) as compared to coalescence results with the DCM model 
(symbols). The small bars on the right hand axis denote results on di- 
baryon yields from a previous RQMD calculation at -^viv = 200 
GeV 1 74] . In addition, the black lines and symbols depict results for 
the production rate of A's from both models, compared to data (grey 
crosses) from (75T.[7^.[77|l. 



Figures [2] and [3] show our results for the mid rapidity 
yields (\y\ < 0.5) of di-baryons and hypernuclei as a 
function of the beam energy Ei„b. In our calculations we 
considered most central (b < 3.4 fm) Pb+Pb/Au+Au 
collisions at Eua> — 1 - 160A GeV. In addition, figure 
|2] shows the A yield (black lines and squares) for the 
two different models compared to data J75J l76l 17711 . In 



these figures, the UrQMD hybrid model calculations 
are shown as lines, while the DCM Coalescence results 
are depicted as symbols. A striking feature of our 
comparison is that, above Ei„b ~ IQA GeV, both 
computations for most (hyper-)nuclei and di-baryons 
agree very well. At lower energies the strange cluster 
production is suppressed in the transport model due 
to the non-equilibrium of strangeness. In the thermal 
calculations restrictions of energy and momentum 
conservation, resulting in a phase space reduction for 
produced strange particles, strongly decreases strange 



particle yields 115 7[ l58l 15911 . This behavior was also 
observed in a core-corona implementation in the hybrid 
model [79]. 



An instructive result is that the yields of most hyper- 
nuclei have a maximum (or saturation) around 10-20 
A GeV of beam energy. Therefore, the investigation of 
hypernuclei can be effectively pursued at these energies. 
On the other hand, the dependence of their yields up to 
energies of ~200 A GeV can help to clarify the mecha- 
nisms of hypernuclei production. 

Noticeably the yields for di-baryons inlcuding S 
hyperons differ strongly with respect to the model 
applied, for the double S state the difference is as 
large as one order of magnitude. The reason for this 
discrepancy can be understood considering that the 
DCM model produces considerably, by a factor of 
5 times, less S's than the UrQMD hybrid model, 
therefore also the dibaryon formation is strongly 
suppressed (note that the experimental S yield is quite 



well reproduced by the UrQMD-hybrid model [80J|79J)- 



Di-baryon production rates have also been calculated 
in a coalescence approach using the RQMD model for 
yTJviv = 200 GeV collisions of Au nuclei B74I1 . To re- 
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Figure 3: Yields per event of different (hyper-)nuclei in the mid ra- 
pidity region (|y| < 0.5) of most central collisions of Pb+Pb/Au+Au. 
Shown are the results from the thermal production in the UrQMD hy- 
brid model (lines) as compared to coalescence results with the DCM 
model (symbols). 



Figure 4: Full acceptance yields per event of different di-baryons cre- 
ated in most central collisions of Pb+Pb/Au+Au. Shown are the re- 
sults from the thermal production in the UrQMD hybrid model (lines) 
as compared to coalescence results with the DCM model (symbols). 



late our calculations to these results, they are indicated 
as the colored bars on the right axis of figure [2] The 
RQMD model used was in particular tuned to reproduce 
multi strange particle yields (such as the S) and the re- 
sults are therefore close to the ones obtained with our 
thermal/hydrodynamic approach. 

Figures [4] and show the integrated {An) yields for 
all considered clusters as a function of beam energy. As 
with the midrapidity results there is a remarkable agree- 
ment between both approaches. However, the integrated 
yields of non-strange nuclei at high energies are system- 
atically larger in the coalescence approach, although the 
mid-rapidity yield was smaller. This observation can be 
explained when the rapidity distribution of the nuclei is 
considered. In the coalescence approach the probability 
to produce a nucleus increases with rapidity and in par- 
ticular in the fragmentation region, where the nucleons 
have small relative transverse momenta and can easily 
coalesce. 

In addition we point out that the coalescence results 
depend on the parameters of the model. As mentioned, 
in the presented results the parameter pc for A's was 
taken equal to the one of the nucleon's. However, the 
hyperon-hyperon and hyperon-nucleon interactions are 



not very well known and we expect that these parame- 
ters may be different for clusters containing A's or even 
E's. In table|2]we demonstrate how the yields of strange 
dibaryon nuclei depend on the momentum parameter 
pc- As discussed previously, we have accordingly re- 
stricted the re parameter, however, by imposing an em- 
pirical limitation related to the nuclear force properties 
that re can not be larger than 4 fm. One can see, we 
expect a very large variation of the yields depending on 
the parameters. For instance, the probability of a bound 
A-nucleon state may decrease by many orders, if we as- 
sume a small pc corresponding to a low binding energy 
of this state. Usually the parameters are fixed by com- 
parison with experiment. Nevertheless, ratios of hyper- 
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Table 2: Dependence of yield of strange dibaryons (per one event) on 
momentum coalescence parameter (pc in units of [MeV/c]), in central 
(b < 3.5/m) Au+Au collisions at 20A GeV 
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Figure 5: Full acceptance yields per event of different (hyper-)nuclei 
created in most central collisions of Pb+Pb/Au+Au. Shown are the re- 
sults from the thermal production in the UrQMD hybrid model (lines) 
as compared to coalescence results with the DCM model (symbols). 



Figure 6: Yields of anti-particle clusters in the mid rapidity region 
(|y| < 0.5) of most central collisions of Pb+Pb/Au+Au as a function 
of y sp/N ■ Shown are only the results from the thermal production in 
the UrQMD hybrid model (lines with symbols). 



nuclei yields should not be changed in the coalescence 
model. 

When the beam energy of the collisions is increased, 
the system created becomes almost net-baryon free. 
This means that the probability to create an anti-particle 
cluster approaches that of the particle cluster. Figure [6] 
shows the results for anti-particle cluster production at 
mid-rapidity fly] < 0.5) in collisions of Pb+Pb/Au+Au 
at center of mass energies of -\JTnn - 3 - 200 GeV. 
We show only results for the UrQMD hybrid model be- 
cause the DCM calculations are restricted to energies 
up to Eiai, - 160A GeV where the statistics needed 
for a meaningful estimate are quite significant. The 
yields of the anti-particle clusters show a monotonous 
increase with beam energy. They show that, at the high- 
est RHIC energy (and at the LHC) the reconstruction of 
^ He might be a feasible task. 

4.1. A special ratio 

In the following we will discuss the double ratio Rh 
defined as: 



R H = A H/ 3 He ■ p/A 



(1) 



for collisions of Pb+Pb/Au+Au and a wide range of 
beam energies. This ratio is especially interesting, as in 



thermal production, it does not depend on the chemical 
potential of the particles (as fugacities cancel), and any 
canonical correction factors for strangeness are can- 
celed. It has been proposed that this ratio is sensitive to 
the local correlation of strangeness and baryon number, 
therefore being a measure of cgs [81]. 



cbs = -3 



(N B Ns)-(N B )(N s ) 
{N 2 S )-(N S ) 2 



(2) 



To calculate Rh we use the above obtained yields for 
hypernuclei and the proton and A yields from the same 
model. For the hadrons the feed down from resonances 
is taken into account as well as the feed down to the 
3 He from the hypertriton. 

Our results for Rh are shown in figure [7] as an excita- 
tion function of the beam energy y/s^N- Rh is evaluated 
for the mid rapidity region of most central (b < 3.4 fm) 
heavy ion collisions. The lines depict results from the 
UrQMD-hybrid model and the symbols denote DCM 
coalescence results. Experimental data are depicted as 
green symbols with error bars. Because experiments 
usually cannot distinguish between A's and E°'s, we 
show Rh in the cases where the A yield includes IP 
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Figure 7: The Strangeness Population Factor Rh = (■ \H/ 3 He) ■ (p/A) 
as a function of ysjvjv for most central collisions of Pb+Pb/Au+Au. 
We compare results from the thermal production in the UrQMD hy- 
brid model (lines) with coalescence results with the DCM model 
(symbols). The red line and symbols denote values of Rfj where the 
A yield has been corrected for the I. contribution. 



Figure 8: Single ratios of(,H/He) (black solid line and circles) and 
\/p (blue dashed line and squares) from the UrQMD hybrid model 
(lines) and DCM model (symbols). 



5. Conclusion 



(black solid line and squares) and where the yield 
is corrected for the S° (red dashed line and circles). 
This is in fact important as there is no experimental 
indication for a bound -LH hypernucleus. 

The double ratio Rh from the hybrid model turns out 
to be almost energy independent. The same behavior 
has been observed in previous thermal calculations 082II . 
On the other hand, the coalescence result increases with 
decreasing beam energy and is in general larger than the 
thermal result. 

To understand this behavior we plotted the single ra- 
tios 3 A H/ 3 He and A/p from our two approaches (lines 
hybrid model and symbols DCM coalescence) in figure 
[8] Here it is obvious that even though the DCM calcu- 
lation produces less A's per proton, the hypernuclei to 
nuclei ratio is still larger. Hence, the A is more likely 
to form a hypernucleus. There seems to be a stronger 
correlation in the transport calculation as in the hydro- 
dynamic description. In fact the qualitative behavior of 
Rh closely resembles the behavior that is expected for 
cbs, the baryon-strangeness correlation, for a hadronic 
gas B83I1 . This observation leads to the conclusion that 
the information on correlations of baryon number and 
strangeness is lost in the thermal calculation because 
here Rh essentially only depends on the temperature. 
On the other hand, in the microscopic treatment the cor- 
relation information survives and Rh captures the trend 
of c BS . 



We have presented results on hyper-nuclei, anti- 
nuclei and di-baryon production in heavy ion collisions 
over a wide beam energy range. To explore the theoret- 
ical uncertainties we applied two distinct approaches: 
firstly, the thermal production with the UrQMD-hydro 
hybrid model and secondly, the coalescence calculation 
within the Dubna hadron cascade model. Concerning 
most hyper-nuclei and di-baryons both approaches 
agree well in their predictions which gives us confi- 
dence in robustness and significance of the obtained 
results. We find that both the non-equilibrium and 
thermal models may be considered as appropriate 
approaches to describe strange cluster production. 
In agreement with previous studies we demonstrate 
that the most promising energy range to produce 
hyper-clusters will be provided by the FAIR and NICA 
facilities, Eyb ~ 10 - 2QA GeV. Anti-matter clusters 
heavier than t are only feasible at RHIC and LHC 
energies. 

The most interesting result of our study is the appar- 
ent difference in the double ratio Rh when we compare 
our thermal results with the coalescence. This differ- 
ence indicates that the information on correlations of 
baryon number and strangeness are visible in the micro- 
scopic coalescence approach, while they are washed out 
in the thermal picture. This could open the opportunity 
to directly measure the strangeness-baryon correlation, 
which may be sensitive to the onset of deconfinement. 
The present status of the experimental data does unfor- 



tunately not allow for a comprehensive comparison with 
our model calculations. We hope that this situation will 
improve in the upcoming RHIC energy scan and FAIR 
experiments. 
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